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Abstract 

We prove that the smooth solutions to the Cauchy problem for the Navier-Stokes 
equations with conserved total mass, finite total energy and finite momentum of 
inertia lose the initial smoothness within a finite time in the case of space of dimen- 
sion 3 or greater even if the initial data are not compactly supported. The cases of 
isentropic and incompressible fluids are also considered. 
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1 System, known results and main problem 



The motion of compressible viscous, heat-conductive, Newtonian polytropic 
fluid in M X M", n > 1, is governed by the compressible Navier-Stokes (NS) 
equations 

dtp+dw,{pu) =0, (1.1) 
dt{pu) + div^(pu 0u) + V^p = DivT, (1.2) 

9t + Pe) + div^ ((^pI^P + + p)u^ = dw{Tu) + kAJ, (1.3) 

where p, u = {ui, ...,Un), p, e, 9 denote the density, velocity, pressure, internal 
energy and absolute temperature, respectively, T is the stress tensor given by 
the Newton law 

T = Tij = p {diUj + djUi) + A divu 5ij, (1.4) 
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where the constants /i and A are the coefficient of viscosity and the second 
coefficient of viscosity, A; > is the coefficient of heat conduction. We denote 
Div and div the divergency of tensor and vector, respectively. We assume that 
/X > 0, A + |/i > 0. 

The state equations liave the forms 

p = Rp6, e = c6, p = Aexp f —] p'^ . (1.5) 



Here A > is a constant, R is the universal gas constant, S = loge — (7 — 
1) logp is the specific entropy, c = ^^j, 7 > 1 is the specific heat ratio. 

The state equations (1.5) imply 

P=(7-l)pe, (1.6) 

which allows us to consider (NS) as a system for the unknown p, u, p. Indeed, 
from (NS) and (1.6) it follows that 

dtp + (u, V^p) + 7pdivu = (7 - 1) E Ti,d,u, + -| A^. (1.7) 

Thus, therefore further we shall consider the system (1.1, 1.2, 1.7), denoted 
(NS*) for short. 

(NS*) is supplemented with the initial data 

{p,u,p)l^ = {po{x),Uo{x),po{x)) e H"^{W), m > [n/2] + 2. (1.8) 



We also consider the isentropic case where the fluid obeys equations (1.1), 
(1.2) and p = Ap'^, we call this system (NSI) for short. 

In the absence of vacuum, the local existence of classical solutions is known. 
Namely, in it is proved that there exist classical solutions, having the Holder 
continuous second derivatives with respect to space variables and the first ones 
with respect to time. In [2] the system of equations of viscous compressible 
fluid is considered as a particular case of combined systems of differential 
equations. The consideration is performed in the Sobolev spaces H"^ with a 
sufficiently large m. The uniqueness of the solution was proved earlier in [3J. 
The existence and uniqueness of local strong solutions in the case where the 
initial density need not be positive and may vanish in an open set were proved 
recently in[l]. 

At the same time there exists a major open problem: to prove or disprove 
that a smooth solution to (NS) (in higher space dimensions) exists globally 



2 



in time. There are partial results concerning the Cauchy problem for (NS) 
away from a vacuum. In [5] it is proved that if there exists a constant p = 
const > such that (po - P, Mo, 5'o) G n L^{W), and the norm ||po - 

p,Uo, SqWh^, fn > [n/^] + 2, is suitably small, then global solution to (NS) 
from C^([0, oo), if™(M"')) exists. In [6J it is shown that for n = 1 the global 
existence takes place without assumptions on the smallness of the norm ||po — 
P, Uq, S'oIIh™-- 

However if the initial density po is compact, then in arbitrary space dimensions 
no solution to (NS) from C^QO, cx)), i/™(M")) exists ([7]). This blowup result 
depends crucially on the assumption about compactness of support of the 
initial density. It does not seem to solve in a negative way the question of 
regularity for (NS). Indeed, (NS) is a model of non-dilute fluids where the 
density is bounded below away from zero, and therefore it is natural to expect 
the problem to be ill-posed when vacuum regions are present at the initial 
time. At the same time, the conservation of mass in the whole space requires 
a decrease of the density down to zero. 

In [7] the author notes that the global smooth solution to (NS) seems to exist 
at least for small data in the case where initial vacuum appears only at infinity. 

Nevertheless, in [8] a sufficient condition for the blow-up in case that the 
initial density is positive but has a decay at infinity was found. In this work 
it was expected a specific time decay of the velocity component that seems 
to be reasonable for a density away from zero. Further, in [9j it was proved 
that if the solution with finite moment of inertia to the (non-heat-conductive) 
Navier-Stokes system is smooth globally in time, then the solution components 
grow as t — s> oo at least as a certain function specific for every component 
4>{t,x), where (p stands either for the density, pressure, velocity or gradient 
of velocity. Namely, we observe all trajectories of particles x{t) that leave in 
a finite time a ball of finite radius Rq and find that provided the solution 
is smooth, the inequality |0(t, a;(t))| < M(t)\x{t)\^, t ^ oo, holds, with a 
continuous function M{t) and a constant A, specific for every (p. Both in [8] 
and [1] the answer whether the solution blows up or not does not depend on 
the initial data, but only on the prescribed decay at infinity. 



Thus, the question remains: is it true that in the case where the support of all 
components of initial data coincides with the whole space the global smooth 
solution exists for any smooth initial data? 

Below we are going to show that, generally speaking, the answer is negative. 
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2 Integral functionals, solution with decreasing components and 
the statement of main theorem 



The system (NS) is the differential form of conservation laws for the material 
volume fl{t); it expresses conservation of mass 

m = j p dx, 
n{t) 

and balance of momentum 

P — j pu dx, 

and total energy 

S= J Qp|M|' + pe) dx ^Ek{t) + Ei{t). 
a{t) 

Here Ek{t) and Ei{t) are the kinetic and internal components of energy, re- 
spectively. 

If we regard Q{t) = M", the conservation of mass, momentum and energy takes 
place provided the components of the solution decrease at infinity sufficiently 
quickly. 



Let us introduce the functionals 

G{t) = ^J p{t, x)\xfdx, F(i) = {u, x)p dx, 

where the first one is the momentum of inertia, the scalar product of vectors 
is denoted as (., .) . 

For technical reasons we impose the following conditions of decay on the so- 
lution components to (NS*) as |a;| — > oo at every fixed t e R+: 

''=o(^). P = .>0. (2.1) 

l«l = ° (^) . = » (^) ■ (2.2) 



4 



If /c 7^ 0, we require additionally 



\De\=ol-^—), |x| ^ oo, teR+. (2.3) 



One can easily verify that these requirements guarantee the conservation of 
the mass m, energy S, momentum P on solutions to (NS*) and ensure a 
convergence of the momentum of mass G{t). 

We impose no restriction on the solution support, however the decay of the 
solution as |x| — i> cx) in the class considered is greater than it is necessary for 
belonging to C\[0, T), if™(M")). 



Definition 1 We will say that a solution {p, u,p) to the Cauchy problem (1.1, 
1.2, 1.1), (1.8) belongs to the class K if it has the following properties for all 
t>0 : 

(i) the solution is classical; 

(a) the solution decays at infinity according to (2.1 - 2.3); 
(Hi) pit, x) > 0; 

(iv) — ^- = o-{t,x) > 0, ||cr(t, x)||^oo(]8n) = o(t°), t ^ oo, where a = 
(^-i)„2+„_2 ,j < 1 + 2 ^ ^ M otherwise. 

Remark 2.1 In particular, {iv) results that S{t,x) > So = inf S{0,x). For 
(NSI) condition {iv) holds trivially. 



Remark 2.2 From (1.7) we have 

n " 

% {dtS + {u, V,S)) = T,,d^u, + kM. 

Evidently, if the right-hand side is non-negative, then S{t, x) > Sq for all 
t > 0. The first item on the right-hand side is non-negative. Indeed, according 

n n n 

TijdjUi = Y 9j{TijUi) - Y UidjTij = 

i,j=l i,j=l i,j=l 

n n n 

2p J2 (dj^^y + A(divn)2 + pY. (M^y + E {^jU,){^^Uj). (2.4) 

i,j=l ij^j i>j 
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For A > this implies the nonnegativity of the right-hand side of (2.4), for 
A < (2.4) can be estimated from below by 

n n n 

(2/i + nX) J2 (^i^i)^ + I^Yl i^j'^if + 2/i XI idjUi){diUj); 

i,j=l i^j i>j 

this expression is nonnegative by assumption on the viscosity coefficients. 

Therefore we can guarantee the uniform boundedness of entropy from below in 
the case k = and in the isothermal case 9 = const. 



Theorem 2.1 Let n > 3, 7 > the momentum P ^ 0. If inf S(0,x) > 
—00, then there exists no global in time solution to (NS*) from the class ^. 

Analyzing condition {iv) we can see that the nondecreasing of entropy along 
trajectories seems natural, whereas the upper bound on the growth of entropy 
can be considered as unreasonable. So, we can re-formulate Theorem 2.1 as 
follows: 

Theorem 2.2 Let us assume > 3, 7 > and P ^ 0. Then any solution 
to (NS*) with properties {i, ii, Hi) such that the entropy does not decrease along 
the particles trajectories, blows up in a finite or infinite time. If the solution 
keeps smoothness for all t > 0, then ||5'||ioo(Kn), ||p||Loo(Mn), ||divM||2,oo(iRn) rise 
at least as 0(t""*"^) as t 00 (the constant a > is indicated in condition 
M.) 

For the isentropic case we have the following version of Theorem 2.1. 

Theorem 2.3 Let n > 3, 7 > jf initial data {po{x) , uo{x)) satisfy that 
the momentum P ^ 0, then the solution to (NSI) from the class A cannot exist 
for all t > 0. 



3 Proof of the theorems 



We begin with the following extension of formula obtained in to the viscous 
case. 

Lemma 3.1 For solutions to (1.1, 1.2, 1.1) with properties (z), (ii) 

G\t) = F{t), (3.1) 

6 



(3.2) 



PROOF. The lemma can be proved by direct calculation using the general 
Stokes formula and taking into account the decay assumption (ii). For exam- 
ple, from (1.1) we get 



= J {u, x)pdx, 

where Sr is a sphere of radius R with the center at the origin. This proves 
equality (3.1). 

Then we get two-sided estimates of G{t). 

Lemma 3.2 //7 < 1 + |, then for solutions with properties (i), (ii) the esti- 



^\ ' Se + F(0)i + G'(O) < G{t) < Se + F(0)i + G'(O) (3.3) 
(for (NS*)) and 



—t^ + F{0)t + G{0) < G{t) < E{0)t^ + F{0)t + G{0) (3.4) 

Zifl 

(for (NSI)) hold. 





mates 



If'j>l + -, then we have 



8t^ + F(0)t + G(0) < G{t) < 



n{7 - 1) 
2 



St'^ + F{0)t + G{0) 



(3.5) 



for (NS*) and 



2m 




e + F{Q)t + G'(O) < G{t) < 



n{l - 1) 
2 



i:(0)i2 + F(0)i + G'(0) 



(3.6) 



for (NSI). 



PROOF. First of all (3.2) result 



G"{t)^2Ek{t)+n{^-l)E,(t) 
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= 2£{t) - (2 - n(7 - l))Ei{t) = n{j - l)S{t) + (2 - n{j - l))Ek{t). (3.7) 

In the case of (NS*) the total energy S is constant, therefore (3.3) and (3.5) 
follow from (3.7) after integration if we take into account nonnegativity of 
Ek{t) and Ei{t). Estimate (3.3) one can find in [10] for the zero viscosity. 

In the case of (NSI) the total energy is only non-increasing, since (1.1) and 
(1.2) for a constant entropy result in 



S'{t) <-vj \Du\^dx < 0, (3. 



with some positive constant v. Then to get the two-sided estimates of G{t) 
from (3.7) we can use the non-increasing of total energy £{t) < £^(0), the 
nonnegativity of Ei{t) and the estimate -Efc(t) > that follows from the 
Holder inequality. This gives (3.4) and (3.6). 



The next step is two-sides estimate of Ei{t). 



Lemma 3.3 For solutions from the class ^ to (NS*) for sufficiently large t 
we have 

Ci C2 

for 7 < 1 + ^, and 

< E,{t) < (3.10) 



Qiy-l)n/2 - ' - G{t)' 

for 7 > 1 + - , with constants Ci, C2 (Ci < C2). The same estimate holds also 
for (NSI) provided P 7^ 0. 



PROOF. The lower estimate is due to [TO]. It follows from the inequality 

27 "(7-1) 

\\f\\L^iR";dx) < <^7,«ll/llL"(R";dx) ll/llLi{R";|x|2da:)' 

together with the lower estimate of internal energy 
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Eiit) > / p\t,x)dx. 

7 — 1 J 



The latter inequality follows from the state equations in (1.5). The constant 

7 — 1 
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with 

n(7-l) -27 



(n + 2)7 — n 



The method of the upper estimate of Ei{t) is also similar to [lU]. Namely, let 
us consider the function Q(t) = 4:G(t)S(t)—F'^(t). The Holder inequality gives 
F2 < 4:G{t)Ekit), therefore £{t) = Ek{t) + Ei{t) > Ei{t) + g|| and 

We notice also that Q{t) > provided the pressure does not equal to zero 
identically. Then taking into account (3.1), (3.2) and (3.7) we have 

Q'{t) = AG'{t)S{t) - 2G\t)G"{t) + AG{t)S\t) = 

= 2(2 - n(7 - 1)) G'{t)Ei{t) + 4G{t)S'{t). (3.12) 

Further, one can see from (3.2) that in the (NS*) case G'{t) > beginning 
from a positive to for initial data, whereas for (NSI) equality (3.2) result in 

p2 

G"{t) > — , 
m 

and we can guarantee the positivity of G'{t) for sufficiently large t only for 
P ^ 0. Thus, for 7 < 1 + f (3.11), (3.12) (and (3.8) for (NSI)) resuh 

Q'(t) ^ 2 - n(7 - 1) G'it) 

Q{t) - 2 Git) ■ ^^-'^^ 

Then (3.13) and (3.11) give 

^'^^)^G(7-iW2(t)' ^ • (3.14) 

If 7 > 1 + ^, then from (3.12) taking into account the lower estimate of Ei{t) 
we get 

Q-(t) < -2((7 - l)n - 2) ^^J;^/,^^y 

Qit) < C' + 4CiG'^-(^-^)"/2^t), C = g(0) -4CiG'i-(^-^)"/2^0), 
and, at last, from (3.11) 
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Thus, since the leading term in the right hand side of the latter inequality is 
the first one, beginning from a moment t > we get 

The proof is over. 



The next lemma is a key point of the theorem's proof. 

Lemma 3.4 Let n > 3, 'j > If \P\ ^ 0, then there exists a positive 

constant K such that for the solutions of the class ^ the following inequality 
holds: ^ 

J \Du\^dx> KE^^^it). (3.15) 



PROOF. First of all we use the Holder inequality to get 



n + 2 
2n 



n-2 
2n 



pudx < / p"+2 dx 



ImI "-2 dx < 



< g 2n(7-l) 



n+2 n-2 

, 2n / \ 2n 

2_ So / f n-2 S 2n \ If, 2n 

e(n + 2)(7-l) cp„ + 2 I I / \u\"-'2 dx 



(3.16) 



Further, using the Jensen inequality we have for — ^ (o^ 7 — 



(7-l)(n + 2) p ^ -v J 

1 r n-2 S 2n \ TOn {'y — 1) EAt) 

_ / e("+2){7-i) c p;7+2 dx < = — L^iZ. 

m J / m mA 



Thus, the latter inequality and (3.16) give 

\P\ <Ki {Ei{t))^^^^ I / \u\'^dx 



n-2 
2n 



(3.17) 



with the positive constant Ki that depends on 'j,n,m,So. Further, we take 
into account of the inequality 



\u\"-'2dx\ < K2 \Du\ dx, 



(3.18) 



where the constant K2 > depends on n, n > 3. The latter inequality holds 
for u G i7^(M") ([H], p.22) and follows from the Sobolev embedding. 
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Thus, the lemma statement, the inequahty (3.15), follows from (3.17), (3.18), 

|p|2 

with the constant K — L . 



We begin from the isentropic case (NSI). 

Proof of theorem 2.3. Prom (3.8), (3.15) we have 

£'{t) <-vK (Ei(t))~^^ . (3.19) 
Together with (3.4), (3.6) and Lemma 3.3 inequality (3.19) implies 

E'{t) < -vKCf^ G'^it) <-Lt'^, 



with a positive constant L, for 7 < 1 + |, and 

. n-2 2(n-2) 

8\t) < -V K G^^(^ {t) < -Li "(7-1), 

for 7 > 1 + ^. In both cases this contradicts to the non-negativity of E{t). 
Thus, theorem 2.3 is proved. 

Proof of theorem 2.1. Let us remind that in the case of (NS*) the total energy 
£ is constant for the solutions of the class ^. However, the derivatives of both 
kinetic and internal components of the total energy can be estimated. Namely, 
taking into account the state equation in (1.5) we have 



dt 

A 



j {u, DivT) dx — j {u, Vp) dx < 



< -1/ y \Du\'^dx J e^ {u, VS) dx - A J e c {u, Vp^) dx, (3.20) 

with a positive constant u. Further, together with (1.1) we obtain 

dEi{t) d ( f Ae-c p^ 



dt dt \J 7 — 1 



dx 



= —. r / e <^p^dtSdx-\ ^ / e<^ p^ ^dtpdx = 

c(7 — 1) J 7 — 1 7 

R" R" 

^(7^ / P'ft - 7^ / P-Miv(p«) dx- 
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/ p^(u,VS)dx < (3.21) 
c(7 - 1) J 

c 7 — 1 y 



c(7 - 1) 



Ei{t) + - J e-o p^{u,VS)dx + A J ec {u,Vp^)dx. 



At last, from (3.20), (3.21), the condition (iv), and estimates (3.9), (3.10), 
(3.3), (3.5) one can get for 7 < 1 + ^ 

Q = j^{Ek{t) + Ei{t)) < 

<-u I \Du\' dx + "^^''^^^j'^-^^^") E,{t) < 

<-uK cf^ {t) + o{e) G-"^ (t) < 
<-Lt^+o(t---M), 
with some positive constant L. Analogously, for 7 > 1 + |, 

= ^(Efe(t)+E,(t)) < -iyKC^G'^{t)+o{r)G-\t) < -Lt'^+oit"-^) 



If we take a from the condition {iv), we get a contradiction that proves the 
theorem. 

Proof of theorem 2.2. It remains to prove the upper estimates for ||p||Loo(Rn) 
and ||divM||ioo(]gn). It is easy to compute that 

E[{t) ^ - j p d\\udx + {n + \) j \d\\u\^ dx + n j \Du\^ dx. (3.22) 

Let us denote /i(t) = ||divM||L°°(R") and /2(t) = ||p||l°°(]R")- As in the proof of 
Theorem 2.1 taking into account Lemma 3.3 we get from (3.22) 

E'^t)>-{^-l)Ei{t)f^{t)+p J \Du\^dx> 
> -(^-l)E,{t)f,{t) + K{E,{t)y^ > -L,t^f,{t) + L2t'^, 
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> 



7-1 



/i + A 



E',{t) > -l^E,{t)f2{t)+fi f \Vu\^dx > 
u + A J 



with positive constants Li, L2, -L3 where /? = — Ti(7 — 1), for 7 < 1 + f , and 
13 = —2, otherwise. If the growth rate of fi(t) and f2(t) is less then prescribed 
in Theorem 3.1 statement, we get a contradiction. Thus, the theorem is proved. 

Remark 3.1 As follows from (3.22) and Lemmas 3.2, 3.3, 3.4, the require- 
ment of incompressibility div^w = signifies that in conditions of Theorem 
2.2 the solution to (NS*) loses its initial smoothness within a finite time. 



3.1 Long-time behavior of solution and the blow up 



In it was found that there exists no global smooth solution to (NS) such 
that 

" {u,x) 



lim sup 

t— >oo 



1 + \x\ 



< 1. 



(3.23) 



L°°(R") 



One can also derive that there exists no global smooth solution to (NS) such 
that 

|P| 



lim sup 



L°°(IR") 



const 



m 



or 



Condition (3.24) follows from the Holder inequality: 



pu dx 



< m \\u\ 



(3.24) 



lim sup / \\divu\\ (t) dr — n\nt \ < const, to>0. (3.25) 



To prove (3.25) (e.g. for 7 < 1 + ^) we note that (3.22) results in 
E'iit) > - J p divudx > -(7 - 1) Ei{t) \\divu\\^^^-^„^ . 

Further, integrating from any to > gives 

t 

In Ei{t) > -(7 - 1) J ||divM||^^(jg„) (r) dr + lnEi(^o)• 



(3.26) 



to 
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Further, inequalities (3.3), (3.9) give the following estimate for sufficiently 
large t: 

Ei{t) < C2 {St^ + F{0)t + ^(O))-^^^^)"/^^ (3 27) 

Estimate (3.25) follows from (3.26) and (3.27) immediately. 

However, condition (3.23) just as conditions (3.24) and (3.25) does not use the 
fact that the velocity belongs to the space H^{M."'). These conditions can be 
applied for the case of zero coefficients of viscosity, i.e. for the gas dynamics 
equations. Moreover, it is possible to construct global in time exact solutions 
to the gas dynamics equations with the velocity that increases by modulus as 
|x| — > 00 and has the form u = A{t)x, A{t) ~ t^^, t ^ 00 (see e.g. [12], [S] 
for details) such that conditions (3.23), (3.24), (3.25) become equalities. These 

solutions satisfy (NS) as well. As follows from (3.25), for smooth solutions to 

t 

(NS) the function / \\divu\\^^(^^„^ (r) dr > nlnt + const, to > 0, and the 

to 

comparison with with the statement of Theorem 2.2 shows that for solutions 
of class H^(R"') this estimate is very far to be exact. 

The author thanks Profs. S.Albeverio and A.A.Zlotnik for a helpful discussion. 
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